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Abstract — In this paper, we study the resource allocation al- 
gorithm design for multiuser orthogonal frequency division mul- 
tiplexing (OFDM) downlink systems with simultaneous wireless 
information and power transfer. The algorithm design is formu- 
lated as a non-convex optimization problem for maximizing the 
energy efficiency of data transmission (bit/Joule delivered to the 
users). In particular, the problem formulation takes into account 
the minimum required system data rate, heterogeneous minimum 
required power transfers to the users, and the circuit power 
consumption. Subsequently, by exploiting the method of time- 
, sharing and the properties of nonlinear fractional programming, 
> the considered non-convex optimization problem is solved using an 
' efficient iterative resource allocation algorithm. For each iteration, 
' the optimal power allocation and user selection solution are 
; derived based on Lagrange dual decomposition. Simulation results 
• illustrate that the proposed iterative resource allocation algorithm 
' achieves the maximum energy efficiency of the system and reveal 
' how energy efficiency, system capacity, and wireless power transfer 
benefit from the presence of multiple users in the system. 

I. Introduction 

Orthogonal frequency division multiplexing (OFDM) is one 
of the leading candidates for supporting high data rate wireless 
broadband communication systems, as envisioned e.g. in the 
, 3GPP Long Term Evolution Advanced (LTE-A) and IEEE 
' 802. 1 1 a/g/n Wireless Fidelity (Wi-Fi) standards, due to its flex- 
ibility in resource allocation and ability in exploiting multiuser 
diversity. In practice, a wireless communication system is ex- 
' pected to support multiple mobile users and to guarantee quality 
of service. However, because of the limited radio resources and 
harsh wireless channel conditions, some of the mobile users are 
typically switched to idle mode since they cannot be served by 
the system temporarily. Unfortunately, mobile devices are often 
battery driven and energy is dissipated even if they are idle 
which creates bottlenecks in perpetuating the network's lifetime. 

Recently, driven by environmental concerns, green mobile 
communication has received considerable interest from both 
industry and academia |fll-|[3l. A promising approach to enhance 
the energy efficiency (bit-per-Joule) of wireless communication 
systems is to harvest energy from the environment. Solar, 
wind, and geothermal are the major renewable energy sources 
for generating electricity. However, these conventional natural 
energy sources may not be suitable for mobile devices and 
not be available in enclosed/indoor environments. On the other 
hand, wireless power transfer, in which energy is harvested from 
propagating electromagnetic waves (EM) in radio frequency 
(RF), is becoming a new paradigm in energy harvesting since 
it recycles the abundant ambient RF energy ||4|-||9|- Although 
the development of wireless power transfer technology is still 
in its infancy, there have been some preliminary applications 
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of wireless power transfer such as wireless body area net- 
works (WBAN) for biomedical implants |0] and passive radio- 
frequency identification (RFID) systems fSj- Indeed, EM waves 
can carry both information and power/energy simultaneously 
ll6l-||9l. The utilization of this characteristic of EM waves 
imposes many new challenges for wireless communication 
engineers. In fSl and fl\, the fundamental trade-off between 
system capacity and wireless power transfer was studied for flat 
fading and frequency selective channels, respectively. However, 
and ||7] assumed a theoretical receiver, which is able to 
decode information and extract power from the same received 
signal but is not yet available in practice. In ||8] and ||9l, 
the authors proposed different power allocation schemes for 
multiple antenna two-user narrowband systems by separating 
the process of information decoding and energy harvesting 
into two receivers. However, if a multicarrier system with an 
arbitrary number of users is considered, the results in [0 and ||9l 
which are valid for single-carrier transmission, may no longer 
be applicable. Besides, the energy efficiency of wireless infor- 
mation and power nansfer systems is still unknown since the 
power dissipations in RF transmission and electronic circuitries 
have not been taken into account in the literature, e.g. |l6]-||9l. 

In this paper, we address the above issues and focus on 
the resource allocation algorithm design for energy efficient 
communication in multiuser OFDM systems with wireless in- 
formation and power transfer In Section |II] we introduce the 
adopted multiuser OFDM channel model. In Section |III1 we 
formulate the resource allocation algorithm design as a non- 
convex optimization problem, which is solved by an efficient 
iterative resource allocation algorithm in Section |IV] Section |V] 
presents simulation results for the system performance, and in 
Section IVH we conclude with a brief summary of our results. 

II. System Model 
In this section, we introduce the OFDM system model. 

A. Multiuser OFDM Channel Model 

We consider a multiuser OFDM system which consists of a 
transmitter and K mobile users. All transceivers are equipped 
with a single antenna, cf. Figure [T] The total bandwidth of 
the system is B Hertz and there are np subcarriers. Each 
subcarrier has a bandwidth W ~ B/np Hertz. We assume that 
the OFDM signaling is time slotted and the length of each time 
(/scheduling) slot is comparable to the length of the channel 
coherence time; the channel impulse response is assumed to 
be time invariant during each scheduling slot. As a result, the 
downlink channel state information (CSI) can be accurately ob- 
tained by exploiting feedback from users in frequency division 
duplex (FDD) systems and channel reciprocity in time division 
duplex (TDD) systems. At the beginning of each scheduling 
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Fig. 1. Illustration of a multiuser OFDM system with K 
downlink wireless information and power transfer 
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slot, the transmitter computes the resource allocation policy 
based on the available CSI. The downlink received symbol at 
user k G {!,..., A'} on subcarrier i S {!,..., rii?} in a 
scheduling slot is given by 



(1) 



where Xi,k, Pi,k, and Hi k are the transmitted symbol, the 
transmitted power, and the multipath fading coefficient between 
the transmitter and user k on subcarrier i, respectively. We 
assume that the transmitted symbol is zero mean with variance 
fjlXifcp} = 1, Vijfc, where £{■} denotes statistical expecta- 
tion. Ik and gk represent the path loss and shadowing between 
the transmitter and user k, respectively. Zi^k is the additive white 
Gaussian noises (AWGN) on subcarrier i at user k with zero 
mean and variance cr^. 

B. Information Decoding and Energy Harvesting Receiver 

In this paper, we assume that each user has the ability to 
decode the modulated information and to harvest energjOl from 
the received radio signal. However, the signal used for decoding 
of the modulated information cannot be used for harvesting 
energy ||8|. As a result, at any given scheduling slot, a user 
can either decode information when it is active (being served 
by the transmitter) or harvest energy when it is idle, but not 
both concurrently, cf. Figure [T] 

III. Resource Allocation 

In this section, we define the system energy efficienc}0 and 
formulate the corresponding resource allocation problem. 

A. Instantaneous Channel Capacity 

In this subsection, we define the adopted system performance 
measure. Given perfect CSI at the user, the channel capacity in a 
scheduling slot between the transmitter and user k on subcarrier 
i with channel bandwidth W is given by 



i.k 



M/log2 1 + P.,fer,,fe and r,,fe = 



hgklHi^k 



(2) 



where Fi fc is the channel-to-noise ratio (CNR) at user k on 
subcarrier i. On the other hand, we assume that in each 
scheduling slot one user is selectec0 for information transfer and 

'The details of the energy harvesting process are beyond the scope of this 
paper and interested readers may refer to I?] for a detailed desciiption. 

^We note that, in the paper, a normalized energy unit is adopted for resource 
allocation algorithm design, i.e., Joule-per-second. Thus, the terms "power" and 
"energy" are used interchangeable. 

' We note that single user transmission in each scheduling slot is commonly 
used in practical multiuser systems such as Wi-Fi systems. On the other hand, 
the adopted framework can be generalized to the case where the data of different 
users is multiplexed on different subcaiTiers, at the expense of a more involved 
notation. 



served by the transmitter. Concurrently, the remaining K — I 
idle users harvest energy from the radio signal emitted by the 
transmitter. Then, the system capacity in a scheduling slot is 
defined as the total number of bits successfully delivered to the 
selected user (bit-per-second) and is given by 



np K 



u{r,s)^}^^wkSka., 



(3) 



i=i k=i 



where V = {P,,k > 0, Vi, k} and 5 = {sfe = {0, 1}, Vfc} are 
the power allocation and user selection policies, respectively. 
Wfc > 0, Vfc, is a positive constant which allows the transmitter 
to give different priorities to different mobile users and to 
enforce certain notions of fairness. On the other hand, for 
designing an energy efficient resource allocation algorithm, we 
incorporate the total power dissipation of the system in the 
optimization objective function. To this end, we model the 
power dissipation (Joule-per-second) in the system as; 

np- K 

Utp{V, S) = + - ^4) 

i=i k=i 

Where Ph = ^^^^'^ ( E ''^ \H^.J?h9^ (5) 



i=i k=i 



Power harvested by idle users 

and j G {!,..., K}. Pq > is the static circuit power 
dissipation of device electronics of the active transceiver such as 
mixers, filters, digital-to-analog converters, and is independent 
of the actual transmitted power. The total power consumption 
of the A' — 1 idle users is omitted in ^ since it is relatively 
small compared to the power dissipation of the active transceiver 
pair. The middle term on the right hand side of (|4) is the 
power consumption in the power amplifier and e > 1 is a 
constant which accounts for the power inefficiency of the power 
amplifier. For example, if e = 5, then 5 Watts are consumed in 
the power amplifier for every 1 Watt of power radiated in the 
RF which results in a power efficiency of i = i = 20%; power 
that is not converted to a useful signal is dissipated as heat in 
the power amplifier. On the other hand, the minus sign in front 
of Ph in (lU indicates that a part of the power radiated by the 
transmitter can be possibly harvested by the idle users. Here, 
we assume that the ability to harvest energy is heterogeneous. 
Specifically, we define the constant parameter < ?7j < 1 
used in ^ as user j's efficiency in harvesting energy from 
the received radio signal. The term 'nj\Hij\'^ljgj in (|5]) can 
be interpreted as frequency selective power transfer efficiency 
for transferring power from the transmitter to idle user j on 
subcarrier i. We note that although UTpiV, S) can be a negative 
value mathematically since \Hij\'^ G [0, oo), UxpiV ,S) > 
always holds in practical communication systems. In particular, 
since the transmitter is the only RF energy source of the idle 
users, we have 2"=! Ef=i ^^^j.fe^fe > Y.t=iJ2k=i Pi,kSk > 
Er=iEf=i^^fcSfc(Ejyfe^jl^»jf where the second 

inequality is due to the second law of thermodynamics in energy 
flow from physic. 

The energy efficiency of the considered system in a schedul- 
ing slot is defined as the total number of bits successfully 
delivered to the selected users per unit energy consumption (bit- 
per-Joule) which is given by 



UeffiV,S) 



Utp{V,S)' 



(6) 



B. Optimization Problem Formulation 

The optimal power allocation policy, V*, and user selection 
policy, S* , can be obtained by solving 



max Ueff{V,S) 



(7) 



s.t. CI: > (l-sfe)^;-^., Vfc, 



.1 j=i 

i=l fc=l 

C3: ^^£P.,feSfe + Pc <^PG, 
i=i fc=i 



A" 



C4: ^ SfcCa: > Rn 

1=1 fe=i 
C6: Sfe G {0,1}, Vfc, 



C5: ^Sfc < 1, 

k=l 

C7: P,,fe > 0,Vi,k. 



Here, ^'^'^^^ in CI denotes the minimum required power transfer 
for user k if it is idle. Pmax in C2 is the maximum transmit 
power allowance to control the amount of out-of-cell interfer- 
ence. C3 constrains the total power consumption of the system 
to not exceed the maximum power supply from the power grid, 
PpG- C4 ensures a minimum required system data rate i?inin. 
C6 is a combinatorial constraint on the user selection variables. 
C5 and C6 are imposed to guarantee that in each scheduling 
slot at most one user is served by the transmitter C7 is the 
non-negative constraint on the power allocation variables. 

IV. Solution of the Optimization Problem 

The optimization problem in (|2ll is a mixed non-convex and 
combinatorial optimization problem. The non-convexity comes 
from the objective function which is the ratio of two functions 
and the combinatorial nature comes from the integer constraint 
for user selection. The first step in solving the considered 
problem is to simplify the objective function using techniques 
from nonlinear fractional programming. 

A. Transformation of the Objective Function 

For the sake of presentation simplicity, we define F as the 
set of feasible solutions of the optimization problem in d?) and 
{V^S} S As the power allocation variables are constrained 
by C2, C3, and C7, is a compact set. Without loss of 
generality, we define q* as the maximum energy efficiency of 
the considered system which is given by 

U{V\S*) U{P,S) 



Utp{P*,S*) v.s Utp{P,S)' 



(8) 



We are now ready to introduce the following Theorem which 
is borrowed from nonlinear fractional programming [flOl . 

Theorem 1: The maximum energy efficiency q* is achieved 
if and only if 



max U{P,S)-q*UTp{P,S) 



v.s 



= U{P*,S*) -q*UTp{P*,S*) ^0, 



(9) 



for U{P,S) > and Utp{P,S) > 0. 

Proof: Please refer to ||3] Appendix A] for a proof of 
Theorem 1. 

Theorem [T] provides a necessary and sufficient condition for 
the optimal resource allocation policy. Specifically, for an opti- 
mization problem with an objective function in fractional form. 



TABLE I 

Iterative Resource Allocation Algorithm. 

Algorithm 1 Iterative Resource Allocation Algorithm 

I: Initialize the maximum number of iterations Lmax and the 

maximum tolerance e 
2: Set maximum energy efficiency q — and iteration index 

71 = 

3: repeat {Main Loop} 

4: Solve the inner loop problem in (flOb for a given q and 

obtain resource allocation policy {P',S'} 
5: if U{P', S') - qUTp{P', S') < e then 
Convergence = true 

UjV'.S') 



return {P* ,S*} = {P' ,S'} q* 
else 



Utp{V',S') 



mv',s') 



and n = n + 1 



1 ^ Utp{V',S') 

Convergence = false 
end if 

until Convergence — true or n = L„ 



there exists an equivalent optimization problem with an objec- 
tive function in subtractive form, e.g. U{P,S) — q*UTp{P,S) 
in the considered case, such that both problem formulations lead 
to the same optimal resource allocation policy. As a result, we 
can focus on the equivalent objective function in the rest of the 
paper. 

B. Iterative Algorithm for Energy Efficiency Maximization 

In this section, an iterative algorithm (known as the Dinkel- 
bach method UOj) is proposed for solving (|7]i with an equivalent 
objective function such that the obtained resource allocation pol- 
icy satisfies the conditions stated in Theorem 1. The proposed 
algorithm is summarized in Table |I] and the convergence to the 
maximum energy efficiency is guaranteed if the inner problem 
([Tol l in each iteration can be solved. 

Proof: Please refer to [IJ] Appendix B] for a proof of 
convergence. 

As shown in Table |I] in each iteration in the main loop, i.e., 
in lines 3-12, we solve the following optimization problem for 
a given parameter q: 



max U{P,S)-qUTp{P,S) 
s.t. CI, C2, C3, C4, C5, C6, C7. 



(10) 



We note that U{P, S)-qUTp{P, S) > holds for any value 
of q generated by Algorithm I. Please refer to 13] Proposition 
3] for a proof. On the other hand, it can be observed that the 
problem formulation in (fTOt for energy efficiency maximization 
is a generalized problem formulation for aggregate weighted 
system capacity maximization. Indeed, if we set q = 0, then the 
objective function in (fTOl i will become the aggregate weighted 
system capacity. 

Solution of the Main Loop Problem: The transformed prob- 
lem is now a mixed convex and combinatorial optimization 
problem. The integer constraint for user selection in C6 is still 
an obstacle in tackling the problem. Indeed, the traditional brute 
force approach or a branch-and-bound method can be used to 
obtain a global optimal solution but result in a prohibitively 
high complexity with respect to (w.rt.) the numbers of users and 
subcarriers. In order to strike a balance between computational 
complexity and optimality, we follow the approach in liTTI and 
relax Sfc in constraint C6 to be a real value between zero and 
one instead of a Boolean, i.e., C6: < Sfe < 1, Vfc. Then, 
Sk can be interpreted as a time-sharing factor for the K users 



to utilize the subcarriers. For facilitating the time-sharing, we 
introduce a new variable and define it as Pi^k = Pi.kSk, Vijk. 
The variable represents the actual transmitted power in the 
RF of the transmitter on subcarrier i for user k under the 
time-sharing assumption. Although the relaxation of the user 
selection constraint will generally lead to a suboptimal solution, 
the authors in |12 | show that the duality gap (suboptimality) due 
to the constraint relaxation becomes zero when the number of 
subcarriers is sufficiently large, e.g. np = 128. 

With this relaxation, it can be shown that the problem is now 
jointly concave w.r.t. the power allocation and user selection 
variables under the time-sharing assumption. As a result, under 
some mild conditions, solving the dual problem is equivalent to 
solving the primal problem 1 1 3 J . 

C. Dual Problem Formulation 

In this subsection, we solve the resource allocation optimiza- 
tion problem by solving its dual for a given value of q. For 
this purpose, we need the Lagrangian function of the primal 
problem in (fTOl i which is given by 

/:(a,/3,7,A,<5,P,5) (11) 

K rip- 
k=l i=l 

np K K 

- ^{Y.Y.^P^'^+Pc-Ppg) -<5(^sfc-l^ 

i=l k=l 

~ ^( X! X! ~ ^max ) - iRn 

i=l k=l 



fe=l 



K np K 

^ a,((l - ..)P;t:5^, -Y^Y. P^^^ {v^\H^A'h9k 



k=l 



1=1 J=l 



Here, a has elements ak. A: G {1, . . . , K}, and is the Lagrange 
multiplier vector accounting for the minimum required power 
transfer for idle users in CI. /3 > is the Lagrange multiplier 
corresponding to the maximum transmit power limit in C2. A > 
is the Lagrange multiplier for C3 accounting for the maximum 
power dissipation in the transmitter due to the limited power 
supply from the power grid. 7 > and S > are the Lagrange 
multipliers associated with the minimum data rate requirement 
and the user selection constraint in C4 and C5, respectively. 
On the other hand, the boundary constraints C6 and C7 on 
the user selection and power variables will be absorbed into 
the Karush-Kuhn-Tucker (KKT) conditions when deriving the 
optimal resource allocation policy in the following. 
Thus, the dual problem is given by 



min max C(a, B.j, X.S.V.S). 

a,/3,7,A,5>0 V.S 



(12) 



D. Lagrange Dual Decomposition 

By Lagrange dual decomposition, the dual problem can be 
decomposed into two layers: Layer 1 (inner maximization in 
(fTST i) consists of np + 1 subproblems where rip of them have 
identical structure and can be solved in parallel; Layer 2 (outer 
minimization in ( fT2] i) is the master problem. The dual problem 
can be solved by solving the problems in Layer 1 and Layer 
2 iteratively, where in each iteration, the transmitter solves the 
subproblems by using the KKT conditions for a fixed set of 
Lagrange multipliers, and the master problem is solved using 
the gradient method. 



Layer 1 (Subproblem Solution): Using standard optimization 
techniques and the KKT conditions, the closed-form optimal 
power allocation on subcarrier i for user k for a given q is 
obtained as 

W{wk + 7) _ 1 "'"^ 



P*k 



Vz,fe,(13) 



where Q 



i.k 



q\e 



J2Vj9jlj\Hr,j\ 

3^k 



Xe 



(14) 



and \x\ = max{0,x}. The optimal power allocation solution 
in ( fT3] l has the form of multilevel water- filling. In particular, the 
water-level, i.e., is different across different subcarri- 

ers and different users. In fact, the water-level on subcarrier i for 
user k depends not only on the priority of user k via Wk, but also 
on its influence on the other K—1 users via J2j^k 'Hj^jh l-^i.j 
Besides, Lagrange multipliers 7 and aj force the transmitter to 
transmit with a sufficient amount of power to fulfill the system 
data rate requirement i?,ni„ and the minimum power transfer 
requirement for idle user j, respectively. 

On the other hand, in order to obtain the optimal user 
selection, we take the derivative of the subproblem w.rt. Sk, 



dC{a,l3.-i,\.S,V,S) 
dsf 



-^kP^sL^, where Qk 



which yields 

is the marginal benefit achieved by the system by selecting user 
k. From ( fT2] i we obtain Qk — 



W{wk + 7) 




Ti,kPlJH2f 



1 



Thus, the optimal user selection is given by 

1 ifQfeb]>o 

otherwise 



(15) 



(16) 



We note that although Pc does not appear in (fT3]l-(fT6]l, it has an 
influence on the solution of the dual problem via the updating 
process of q, cf. Table I. 

Layer 2 (Master Problem Solution): For solving the Layer 2 
master minimization problem in (fT2t . i.e, to find a, /3, 7, A, and 
6 for given V and S, the gradient method can be used since the 
dual function is differentiable. The gradient update equations 
are given by: 



np K 



ak{m + l)=\ak{m) - Ci(m) x {j_^2^Pi 



«J ^3 



(vklHi^kl'^lkgkj 



np K 
i=l fc=l 



7(m + 1)= 7(m) - ^sim) x 
A(m + 1)= A(m)- ^4(771) X 
-Pc 



np K 



2^ 2^ SkC^,k 



Ci.k Ru 



1=1 fc=l 
PpG 



np- K 

EE'^kP^k 

1=1 fc=i 



6{m + 1)= S{m) - ^5 (to) X 1 - 



K 

fc=i 



(17) 
(18) 



t(19) 



(20) 



(21) 
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Fig. 2. Average system energy efficiency (bit-per- Joule) versus number of 
iterations with different numbers of users, K, and different values of maximum 
transmit power allowance, Pmax- 

where index m > is the iteration index and ^„(m), 
u S {1, . . . , 5}, are positive step sizes. The updated Lagrange 
multipliers in (fT7]i-(l2ni are used for solving the subproblems 
in (fT2] i via updating the resource allocation policy in (fT3Tl-(fT6]l. 
Since the transformed problem is concave for a given parameter 
q, it is guaranteed that the iteration between Layer 2 (master 
problem) and Layer 1 (subproblems) converges to the primal 
optimum of (fTOl i in the main loop, if the chosen step sizes 
satisfy the infinite travel condition ifTSl . 

To summarize the iterative algorithm between Layer 1 and 
Layer 2, the gradient update in (fT7i- (l2]l can be interpreted as 
the pricing adjustment rule of the demand and supply model 
ifTJl . Specifically, the Lagrange multipliers can be interpreted 
as a set of shadow prices for utilizing the resources. If the 
demand of the resource exceeds the supply, then the gradient 
method will raise the shadow prices via adjusting the Lagrange 
multipliers in the next iteration; otherwise, it will reduce the 
shadow prices until some users can afford them. By combining 
the gradient update equations and the user selection criterion 
in ( fTSI l, only one user is selected eventually even though time- 
sharing is introduced for solving the transformed problem in 
(Ell. 

V. Simulation Results 

In this section, we evaluate the performance of the proposed 
resource allocation algorithm using simulations. An indoor 
communication system with a maximum service distance of 10 
meters is considered. The TGn path loss model [il4< is adopted 
with a reference distance of do = 2 meters. The K users are 
uniformly distributed between the reference distance and the 
maximum service distance. The effective antenna gain for each 
transceiver is assumed to be 14 dB. The system bandwidth is 
S = 5 MHz, the number of subcarriers is np = 128, and 
Wk = 1, Vfc. Note that by setting Wk = 1, Vfc, we obtain 
the maximum achievable system energy efficiency. We assume 
a carrier center frequency of 470 MHz which will be used 
by IEEE 802.11 for the next generation of Wi-Fi systems 
ifTSl . Each subcarrier for RF transmission has a bandwidth of 
W = 39 kHz and the noise variance is <t^ = —118 dBm. 
The multipath fading coefficients of the transmitter-user links 
are generated as independent and identically distributed (i.i.d.) 
Rayleigh random variables with unit variance. The shadowing 
of all communication links is set to dB, i.e., gk = 1, Vfc. 
We assume a static circuit power consumption of Pc = 40 
dBm, a maximum power grid supply of Ppc = 50 dBm, and a 
minimum data rate requirement of i?,nin = 10 Mbits/s. The 
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Fig. 3. Average system energy efficiency (bit-per-Joule) versus maximum 
transmit power allowance, Pmax, for different numbers of users, K. 

minimum required power transfer and the energy harvesting 
efficiency are set to -P^tnt ~ dBm, Vfc, and rjk = 0.8, Vfc, 
respectively. Besides, we assume a power efficiency of 40% for 
the power amplifier used at the transmitter, i.e., £ = jf^ ~ 2.5. 
The average energy efficiency is computed according to (|6) and 
averaged over multipath fading and path loss. In the sequel, the 
total number of iterations is defined as the number of main 
loops in Algorithm 1 multiplied with the number of iterations 
in solving the Layer 1 and Layer 2 problems. Moreover, the 
step sizes adopted in (fT7]i- (l2Tl i are optimized for obtaining a 
fast convergence. Note that if the transmitter is unable to fulfill 
the minimum required system data rate i?min or the minimum 
required power transfer P^^i^^k' energy efficiency 

and the system capacity for that channel realization to zero to 
account for the corresponding failure. 

A. Convergence of Iterative Algorithm 1 

Figure |2] illustrates the convergence behavior and the duality 
gap of the proposed iterative algorithm for maximizing the 
system energy efficiency. The duality gap is defined as the dif- 
ference between primal optimum and dual optimum (achieved 
by the proposed algorithm) and indicates the suboptimality due 
to the constraint relaxation in C6 and insufficient numbers of 
iterations. We investigate the system performance for different 
numbers of users, K, and different values for the maximum 
transmit power allowance, Pmax- The results in Figure |2] were 
averaged over 10^ independent channel realizations for both 
path loss and multipath fading. It can be observed that the 
energy efficiency of the proposed algorithm converges to the 
optimum value within 30 iterations for all considered scenarios. 
The fast convergence and the zero duality gap confirm the 
practicality of the proposed algorithm. 

In the following results, we set the total number of iterations 
to 30 to illustrate the performance of the proposed algorithm. 

B. Average Energy Efficiency and Average System Capacity 
Figure [3] depicts the average system energy efficiency versus 

the maximum transmit power allowance, Pmax, for different 
numbers of users, K. It can be observed that the energy 
efficiency of the proposed algorithm increases w.rt. Pmax 
monotonically and reaches an upper limit where the energy 
efficiency gain due to a higher value of Pmax vanishes. This 
result indicates that once the maximum energy efficiency is 
achieved by transmitting a sufficiently large power in the RF, 
any additional increase in the transmitted power will incur a loss 
in energy efficiency. On the other hand, the energy efficiency 
of the system increases with the number of users and there are 
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Fig. 4. Average system capacity (bit-per-second) versus maximum transmit 
power allowance, Pmax, for ditferent numbers of users, K. 

two reasons for this behaviour First, as the number of users 
in the system increases, the transmitter has a higher chance to 
select a user who has a strong channel due to multiuser diversity 
(MUD). Indeed, the MUD introduces an extra power gain to 
the system which helps save energy. Second, more idle users 
harvest the power radiated by the transmitter which reduces the 
total power consumption of the system, cf. (|4|i. For comparison, 
Figure|3]also contains the energy efficiency of a baseline scheme 
which adopts a resource allocation algorithm maximizing the 
system capacity (bit/s) under constraints C1-C7. It can be 
seen that in the low transmit power allowance regime, the 
proposed algorithm performs virtually the same as the baseline 
scheme. Indeed, the small power radiated by the transmitter 
creates a bottleneck in the system and the performance of the 
proposed algorithm is restricted by the limited system resources. 
However, in the high transmit power allowance regime, the 
energy efficiency of the baseline scheme decreases dramatically 
since an exceedingly large transmit power is used for capacity 
maximization. 

Figure m shows the average system capacity versus maximum 
transmit power allowance, Pmax, for different numbers of users, 
K. The proposed algorithm is compared with the baseline 
scheme described in the last section. It can be observed that 
both schemes benefit from an increasing number of users due 
to MUD. On the other hand, the proposed algorithm achieves 
virtually the same system capacity as the baseline scheme in 
the low transmit power regime. This suggests that the proposed 
algorithm transmits with full power in the low transmit power 
allowance regime. However, as the transmit power allowance 
increases, the baseline scheme outperforms the proposed algo- 
rithm, since the former scheme uses all the available power 
for capacity maximization which impairs the system energy 
efficiency. 

C. Average Harvested Power 

Figure |5] shows the total average power harvested by the idle 
users versus the maximum transmit power allowance, Pmax, for 
different numbers of users, K. It can be observed in Figure |5] 
that the total average harvested power increases with Pmax since 
more power is available in the RF. On the other hand, a larger 
portion of the radiated power can be harvested when there are 
more users in the system since more idle users are participating 
in the energy harvesting process. 

VI. Conclusions 
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Fig. 5. Average total harvested power versus maximum transmit power 
allowance, Pmax, for different numbers of users, K. The double-sided an'ows 
indicate the power harvesting gain due to an increasing number of users. 

and combinatorial optimization problem, in which concurrent 
wireless information and power transfer were considered. By 
exploiting nonlinear fractional programming and Lagrange dual 
decomposition, a novel iterative resource allocation algorithm 
was proposed for maximizing the system energy efficiency. Sim- 
ulation results showed that the proposed algorithm converges 
within a small number of iterations and unveiled the potential 
benefits of having multiple users for energy efficiency, system 
capacity, and wireless power transfer. 
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In this paper, we formulated the resource allocation algo- 
rithm design for multiuser OFDM systems as a non-convex 



